Signatures of non-monotonic d— wave gap in electron-doped cuprates 
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We address the issue whether the data on optical conductivity and Raman scattering in electron- 
doped cuprates below T c support the idea that the d— wave gap in these materials is non-monotonic 
along the Fermi surface. We calculate the conductivity and Raman intensity for elastic scatter- 
ing, and find that a non-monotonic gap gives rise to several specific features in optical and Ra- 
man response functions. We argue that all these features are present in the experimental data on 
Nd2- a; Ce; C Cu04 and P^-xCe^CuC^ compounds. 
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I. INTRODUCTION 

The studies of electron-doped cuprates, 
Nd 2 _ 2; Ce a; Cu04- I 5 (NCCO) and Vv 2 - x Ge x CnO A ^ s 
(PCCO) are attracting considerable attention from high- 
T c community. The phase diagram of electron-doped 
cuprates is not as involved as in hole-doped materials. 
It contains sizable regions of antiferromagnetic and su- 
perconducting phases, and only a small region showing 
pseudogap behavior—. The superconducting dome is 
centered around a quantum-critical point at which the 
antiferromagnetic Tn vanishes, in close similarity to 
phase diagrams of several heavy- fermion materials 2 . 

Scanning SQUID 3 and ARPES experiments^ on 
electron-doped cuprates provided strong evidence that 
the gap symmetry is d x 2_ y 2, same as in hole-doped 
cuprates. This gap has nodes along the diagonals of the 
Brillouin zone, and changes sign twice along the Fermi 
surface. The functional form of the d x i_ y i gap is a more 
subtle issue, however. In hole-doped cuprates, the gap 
measured by ARPES follows reasonably well a simple 
d— wave form A(fc) = ^-(cosfe x — cosk y ) (equivalent to 
cos 20 for a circular Fermi surface), at least near and 
above optimal doping^. In the electron-doped cuprates, 
high-resolution ARPES data on the leading-edge gap in 
Pr .89LaCeo.iiCuC>4 (Ref^) show a non-monotonic gap, 
with a maximum in between nodal and antinodal points 
on the Fermi surface. Such gap was earlier proposed in 
Refi^ as a way to explain Raman experiments in NCCO, 
particularly the higher frequency of the pair-breaking, 
'2A' peak in the B 2g channel than in the B\ g channel. 
Recent measurements of optical conductivity <Jx(uS) in 
Pri. 85 Ceo.i5Cu04 (Ref£) were also interpreted as an in- 
direct evidence of a non-monotonic gap. 

The interpretation of the experimental results is still 
controversial, though. ARPES data on PCCO below T c 
in Ref4 show a non-monotonic leading-edge gap, but the 
spectral function all along the Fermi surface does not dis- 
play a quasiparticlc peak, from which one would generally 
infer the functional form of the gap more accurately. The 
interpretation of the Raman data has been criticized in 



RefJi on the basis that, within BCS theory, the shapes 
of B\g and B 2g Raman intensities for the non-monotonic 
gap proposed in Ref. do not agree with the data. Fi- 
nally, optical results for PCCO in Ref. [H do show a max- 
imum at about 70meV, which is close to 2A max inferred 
from Big Raman scattering. However, it is a-priori un- 
clear whether one should actually expect such maximum 
in the optical conductivity. In particular, in hole-doped 
materials, cti(w) is rather smooth at 2A maJ r^. 

From theory perspective, the non-monotonic d— wave 
gap appears naturally under the assumption that the 
d x 2_ y 2 pairing is caused by the interaction with the 
continuum of overdamped antiferromagnetic spin fluc- 
tuations. Spin-mediated interaction is attractive in the 
d x 2_ y 2 channel and yields a gap which is maximal near 
the hot spots - the points along the Fermi surface, sep- 
arated by antiferromagnetic momentum, Qa_f- In opti- 
mally doped NCCO and PCCO, hot spots are located 
close to Brillouin zone diagonals, and one should gener- 
ally expect the d x 2_ y 2 gap to be non-monotonioii. More 
specifically, in the spin fluctuation scenario, the maxi- 
mum of the gap is slightly shifted away from a hot spot 
towards antinodal region, such that d— wave supercon- 
ductivity with a non-monotonic d x 2_ y 2 gap survives even 
when the hot spots merge at the zone diagonals 12 . The 
solution of the gap equation in this case yields 

A - A ( 2 ^ \ cos2< t> m 

mra Uvay (i + acos 2 2,/>) 3 / 2 - w 

Here, is the angle along the (circular) Fermi surface 
(0 = 7r/4 corresponds to a diagonal Fermi point), a > 1/2 
is a model-dependent parameter, and A max is the maxi- 
mum value of the gap located at cos 2(/> = (1/20) 1 / 2 . The 
gap at various a is shown in Fig[TJ As a increases, the 
nodal velocity increases, the maximum of the gap shifts 
towards the zone diagonal, and the value of the gap at 
the antinodal point = decreases. A similar functional 
form of the gap can be obtained by adding higher har- 
monics cos 60, cos 100 etc. to the cos 20 gap. We have 
found, however, that Eq. |T]) is somewhat better for ex- 
perimental comparisons than the gap with a few higher 
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FIG. 1: (color online) The non-monotonic d ;c 2_ a 2-wave gap 
calculated using Eq.{l]) for a = 2 (dotted curve), a — 4 
(dashed curve), and a — 20 (dashed-dotted curve). The 
black straight curve is a monotonic d x 2_ y 2-wa,ve gap, A^> = 
A max cos 20. 



harmonics. ARPES measurements^ place the maximum 
of the gap slightly below <f> = ir/6. This is best repro- 
duced if we set a = 2. However, since ARPES results 
have not been yet confirmed by other groups, we will 
keep a as a parameter and present the results for various 
values of a. 

The goal of our work is to verify to which extent op- 
tical conductivity cr 1 (w) and Raman scattering R(lo) in 
a d— wave superconductor with a gap given by Eq. (TT|) 
are consistent with the experimental data. For this, we 
computed <j\ (u) and R(u>) in B lg and B 2g geometries as- 
suming that the scattering is elastic. The latter does not 
necessarily have to come from impurities - scattering by 
collective excitations in spin or charge channels is also 
dominated by processes with small frequency transfers. 
For simplicity we assume that the normal state damping 
rate is independent on frequency and only focus on the 
effects associated with the pairing. 

We have found several features which distinguish opti- 
cal and Raman responses in superconductors with non- 
monotonic G?-wave gap from superconductors with a 
cos 20 gap. Optical conductivity in a pure d— wave super- 
conductor with elastic scattering has a weak maximum 
followed by a broad suppression region at frequencies of 
the order A max (RefJ^). For the non-monotonic gap, we 
have found a rather strong maximum in ai(u>) slightly 
below 2A max , followed by a sharp drop in conductiv- 
ity down to very low frequencies, where the conductiv- 
ity begins to increase again towards a constant value at 
uj = 0+ (see Fig. For Raman scattering, we have 
observed that the peak in the B 2g channel is located at 
a higher frequency than in the B\ g channel, and also 
that the shapes of the two Raman profiles are very dif- 
ferent - the Bi g peak is near-symmetric, the B\ g peak is 
very asymmetric with shoulder-like behavior above the 
peak frequency. We argue that these features are con- 



sistent with the experimental conductivity and Raman 
data. From this perspective, our findings give additional 
support to the idea that the d x ?_ y i gap in electron-doped 
cuprates is highly non-monotonic. 

We present the formalism in Sec. [Til an d the results 
in Sec. IIIII In the latter section, we also consider the 
comparison with the data in more detail. The last section 
is the conclusion. 



II. THE FORMALISM 

We adopt a conventional strategy of analyzing opti- 
cal and Raman responses in non-s— wave superconduc- 
tors with impurity scatterin g 13 i 14 i 15 . We assume that 
the scattering originates from the s— wave component of 
the effective interaction (which includes the impurity po- 
tential), and gives rise to a k— independent fermionic self- 
energy E(w). The pairing comes from a different, d— wave 
component of the interaction. As in earlier work a 13 i 14 ' 15 , 
we assume that the d— wave anomalous vertex is fre- 
quency independent, and to a reasonable accuracy can be 
replaced by A(0) from Eq. (fTJ). The time-ordered nor- 
mal and anomalous fermionic Green's functions in this 
approximation are given by 
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where Q = uj+Y,{uj). The self-energy is by itself expressed 
via the (local) Green's function via 
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where the local Green's function is 
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( Gl = 1 in the normal state), and the parameter C 
is interpolated between C >> 1 in the Born limit, and 
C << 1 in the unitary limit. 

Optical conductivity <7i(w) and Raman intensity R(u>) 
are both given by the combinations of bubbles made out 
of normal (GG) and anomalous (FF) Green's functions. 
Optical conductivity is proportional to the current- 
current correlator, while Raman intensity is proportional 
to the density-density correlator weighted with angle- 
dependent Raman vertex factors 



7s oc cos 20, 7 B2g oc sin 2<p. 



(6) 



To a first approximation, B\ g Raman scattering then 
gives information about electronic states in the antinodal 
regions, near = 0, while Bi g scattering gives informa- 
tion about nodal regions, near <j) — 7r/4. 
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The overall sign of the FF contribution is different 
for <j(lo) and R(lo), as the running momenta in the side 
vertices in the FF term are k and — k, between which the 
current operator changes sign, but the density operator 



remains intact. For a constant density of states, which 
we assume to hold, the integration over in GG and 
FF bubbles can be performed exactly, and yields 
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Here, o\ is the real part of the conductivity, lo v i is the 
plasma frequency, the index i labels the various scatter- 
ing geometries, Co = lo + S(w), lo± = lo ± ~, and Rq is 
the normalization factor for the Raman intensity. The 
conductivity in a superconductor also contains a <$(f2) 
contribution (not shown) related to the superconducting 
order parameter. 

In a ideal BCS superconductor with Co = lo, the con- 
ductivity ai(uo) vanishes, while Raman intensity is given 
bylkil 



Ri{u>) = RqRc i 
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where < ... > denotes the averaging over the Fermi sur- 
face. For a pure d— wave gap, B\ g Raman intensity scales 
as lo 3 at small frequencies^, and diverges logarithmically 
at 2A max . Big intensity scales as lo at small frequencies 
and has a broad maximum at around \.§A max (see Fig. 
[5] and [5] below) . 
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FIG. 3: The behavior of ImS(w) in the unitary limit. Solid 
curve - a monotonic d— wave gap, dashed and dashed-dotted 
curves - a non-monotonic gap from Eq. ([lj with a — 4 and 
a — 2, respectively. We used 7 = 0.3A. 



III. THE RESULTS 
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FIG. 2: The behavior of /mi = Im(uj + E(a>) = ZmS(ii)) in 
the Born limit for a d-wave superconductor with a monotonic 
gap (solid curve) and non-monotonic gap from Eq. ([lj with 
a = 2 (dashed-dotted curve) and a = 4 (dashed curve). We 
used 7 = 7/C 2 = 0.05A. 



A. Fermionic self-energy 

We computed the fermionic self-energy by solving nu- 
merically the self-consistent equation ((4]) in the Born and 
unitary limits. The results for the imaginary part of 
the self-energy in the Born limit are presented in FigfJl 
For a d— wave superconductor with a monotonic gap, 
ImE is linear in frequency at small lo and has a cusp 
at to = A max . For a non-monotonic gap, ImE is re- 
duced at small frequencies, and then rapidly increases 
to a value comparable to that for a monotonic gap. 
This behavior resembles, particularly for a = 2, the for- 
mal solution of Eq. Q for an angle-independent gap 
A = A max . In the latter case, ImE = up to a frequency 
lo = A(l — (7/A) 2 / 3 ) 3 / 2 , where 7 = 7/C 2 , and rapidly 
increases above this frequency. For 7 = 0.05 A max , used 
in Fig. [21 the jump occurs at u> « 0.8A, much like in the 
plot for a = 2. We emphasize that the solution of ([J| for 
a constant A is not the result for an s-wave supercon- 
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ductor. For the latter, the fermionic self-energy and the 
pairing vertex are renormalized by the same interaction, 
and the self-consistent equation for £(w) does not have 
the form of Eq. ^ with frequency independent A. 

In Fig. [3]we show ImS in the unitary limit C = 0. We 
observe the same trend. For a monotonic <i-wave gap, 
ImE is nearly monotonic, and has only a slight minimum 
around Q.8A max . For a non-monotonic gap, particularly 
for a = 2, ImE has a more pronounced structure with 
a sharp minimum around 0.8A max . This behavior again 
resembles that for a constant gap A. In the latter 
formal solution of ([JJ for 7 <C A yields a zero ImS(w) be- 
tween 2y/jA and \/ A 2 + r y 2 . At larger id, ImS gradually 
approaches the normal state value 7, at small frequencies 
it is also finite and approaches y/jA at zero frequency 
[for a generic C, a non-zero imE(cJ = 0) (the unitary 
resonance) appears when 7 exceeds AC 2 Vl + C 2 ]. The 
region of vanishing ImS shrinks to zero when 7 exceeds 
the critical value of 2A/(3\/3) w 0.4 A. For the same 
7 = 0.3A max as used in Fig. [3J ImE for a constant 
gap sharply drops around l.lA max , and rebounds both 
at larger and smaller frequencies, much like our actual 
solution for a = 2. 




FIG. 4: The behavior of the optical conductivity <ri(u) in 
the Born (a) and unitary (b) limit, for a monotonic gap 
A(<f>) oc cos 20 (solid line), and non-monotonic gap with a = 2 
(dashed-dotted line) and a — 4 (dashed line). 



B. Optical conductivity 

Substituting the results for the self-energy into Eq , 
we obtain optical conductivity. The results are plot- 



ted in Fig. [4] (a) and (b) for Born and unitary limits, 
respectively^. The behavior of the conductivity in the 
two limits is not identical, but the interplay between the 
monotonic and the non-monotonic gap is similar. In 
both cases, the conductivity for a non-monotonic gap 
passes through a well-pronounced maximum at some fre- 
quency below 2A rnax , sharply drops at smaller frequen- 
cies, and then increases again at very low frequencies, and 
at to — > + approaches the universal limit 13 in which the 
conductivity depends on the nodal velocity but does not 
depend on 7 as long as 7 << A max . The universal be- 
havior is, however, confined to very low frequencies, while 
in a wide frequency range below 2A max the conductivity 
in case of a non-monotonic gap is strongly reduced com- 
pared to its normal state value. The frequency at which 
the conductivity has a maximum depends on a, and is 
closer to 2A max for a = 2 than for a — 4. 

The existence of the maximum in <7i(w) below 2A max 
can be also understood analytically. Expanding the gap 
A(<f>) near its maximum value A max and substituting 
the expansion into ((7]), we find, after some algebra, that 
the conductivity has a one-sided non-analyticity below 
uj = 2A max - it contains a negative term proportional 
to (2A max — a;) 3 / 2 . This negative term competes with a 
regular part of a\(uj), which smoothly increases with de- 
creasing u>, and gives rise to a maximum in o\ (uj) below 
2A 

The behavior of the conductivity in a superconductor 
with a non-monotonic gap is consistent with the available 
data on ai(u>) in optimally doped PCCO^. The mea- 
sured conductivity has a rather strong peak at 70cto _1 , 
and decreases at smaller frequencies. The authors of Ref. 
[H explained the existence of the maximum in the opti- 
cal conductivity by a conjecture that the conductivity 
in a d— wave superconductor with a non-monotonic gap 
should largely resemble the conductivity in an s— wave 
superconductor. Our results are in full agreement with 
this conjecture. The authors of Ref£ also associated the 
peak frequency with 2A max . We found that the peak 
frequency is actually located below 2A max , and the dif- 
ference between the two depends on the shape of the gap. 
For our a — 2, the peak frequency is at l.8A max in the 
Born limit, and at 1.3A max in the unitary limit. For 
a = 4, the deviations are higher. Experimentally, 2A max 
in optimally doped PCCO can be extracted from Bi g Ra- 
man scattering (see below) and equals 77cto _1 , see Ref. 
[H, i.e., the peak in o\(u)) is at 1.8A max . This agrees 
with our a — 2 case in the Born limit. 



C. Raman intensity 

The results for the Raman intensity are presented in 
Figs. [S][ni In a BCS superconductor with a monotonic 
gap, Raman intensity has a sharp peak at 2A max in B\ g 
scattering geometry, and a broad maximum at around 
\.&A max for Big scattering. This behavior holds in the 
presence of impurity scattering, both in Born and unitary 
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limits, see Figs. [SJa) and[6ja). 

The Big and Bi„ Raman intensities for 



>2g 



monotonic gap are presented in Figs. -H] (b)-(c) for 
Born and unitary limits and a — 2 and a = 4. In all 
cases, we find the opposite behavior: B2 g intensity has a 
sharp peak at 2A max , while B\ g intensity is very small 
at small frequencies, rapidly increases around A max , 
passes through a maximum, then gradually decreases at 
higher frequencies and displays a weak kink-like feature 
at 2A max . The position of the Bi g peak depends on a 
- in both Born and unitary limits it is close to 1.6A max 
for a = 2, and is close to A r „„ r for a = 4. 
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FIG. 5: The behavior of the Raman intensity R{ui) in Bi g 
(dashed) and Big (solid) scattering geometries for a mono- 
tonic gap A(4>) oc cos 20 (a), non-monotonic with a=2 (b) 
and a=4 (c) in the Born limit. 

The occurrence of the '2A' peak in Bi g channel at a 
higher frequency than in the B\ g channel was the main 
motivation in Refii to propose a non-monotonic d— wave 
gap. The argument was that the gap with a maximum at 
intermediate < 4> < 7r /4 will have more weight in the 
nodal region and less in the antinodal region, thus in- 
creasing the effective '2A' for Bi g intensity and decreas- 
ing it for Big intensity. In the optimally doped PCCO, 
the B^ g peak occurs at 77cm _1 , while the maximum in 
Bi g scattering is around 60cm -1 . In optimally doped 
NCCO, the Big peak occurs at 67cm -1 , while the max- 



imum in Big scattering is at 50cm -1 (Ref. 19. The 
ratios of the peak positions are 1.28 in PCCO and 1.34 
in NCCO. This is consistent with our result for a = 2 
(same a that gives the best fit of ARPES and conduc- 
tivity data), for which this ratio is 1.25. Also, tak- 
ing experimental 67cm -1 for 2A in NCCO, we obtain 
A max = A.2meV , in reasonable agreement with 3.7meV 
observed in tunneling 20 . In addition, the data in Fig. 
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FIG. 6: The behavior of the Raman intensity R{ui) in B\ g 
(dashed) and Big (solid) scattering geometries for a mono- 
tonic gap A((j>) oc cos 20 (a), non-monotonic with a=2 (b) 
and a=4 (c) in the unitary limit 

3 of Ref. 7 show that the B^g peak is nearly symmetric, 
while Big intensity is asymmetric around a maximum - 
it rapidly increases at frequencies around 40cm -1 , passes 
through a maximum, and then gradually decreases at 
higher frequencies. This behavior of Bi g intensity is fully 
consistent with Fig. EJ-H] (b)-(c). 

Blumberg et &\£ also analyzed Raman intensities at 
various incident photon frequencies and found resonance 
enhancement of the Bi g intensity, but no resonance en- 
hancement of Big intensity. We didn't attempt to an- 
alyze the resonance behavior of the Raman matrix ele- 
ment (this would require to consider the internal com- 
position of the Raman vertes 2 ^). We note, however, 
that the shape of the Big Raman intensity virtually does 
not change between the resonance and the non-resonance 
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cases, only the overall magnitude increases near the reso- 
nance much like it happens in resonant Raman scattering 
in insulating cuprates^. We therefore believe that our 
analysis of the Raman profile as a function of transferred 
frequency is valid both in the non-resonance and in the 
resonance regimes. 

Finally, we note that our results for R(lu) are quite sim- 
ilar to Refi^, whose authors criticized the explanation of 
the Raman data in terms of non-monotonic gap2i. How- 
ever, contrary to Ref. — , we argue that the theoretical 
results for R(to) obtained for a non-monotonic gap agree 
well with the data in Refi At the same time, we agree 
with Ref. — that one can hardly extract from the data 
in NCCO and PCCO the lu 3 behavior of B lg intensity 
(which is the Raman hallmark of d x 2_ y 2 pairing), as the 
low-frequency behavior of the B\ g intensity is dominated 
by a sharp increase at frequencies of order A max . 

In the analysis above we neglected final state interac- 
tion (the renormalization of the Raman vertex). There 
are two reasons for this. For Bi g scattering, the final 
state interaction is given by the B^g component of the 
effective four-fcrmion interaction. This component is re- 
pulsive, at least if the effective four-fermion interaction 
comes from spin-fluctuation exchange. The repulsive fi- 
nal state interaction does not give rise to excitonic reso- 
nances, and generally does not substantially modify the 
Raman profile 24 . For B\ g scattering, final state inter- 
action is the same as the pairing interaction, i.e. it is 
attractive. In general, such interaction affects the Ra- 
man profile 2 ^. However, the interaction which gives rise 
to a non-monotonic gap in the form of fl} is the largest 
at angles close enough to tt/4. At these angles, the 
Big matrix element 7s lg cx cos 2<f> is reduced, and we do 
not expect that repeated insertions of B\ g vertices will 



substantially modify the Raman profile. 



IV. CONCLUSION 

In this paper, we analyzed the behavior of the opti- 
cal conductivity and Raman intensity in B\ g and Bi g 
scattering geometries in the superconducting state of 
electron-doped cuprates. We found that the results are 
best fitted by a non-monotonic d x 2_ y 2 gap. Such gap was 
originally suggested as a way to explain Raman data£, 
and later extracted from ARPES measurements of the 
leading edge gap along the Fermi surface^. The non- 
monotonic gap has also been obtained theoretically in 
the analysis of quantum-critical pairing mediated by the 
exchange of overdamped spin fluctuations^ 2 -. 

We found that the non-monotonic gap which agrees 
best with the ARPES data (Eq. (P) with a = 2) also 
fits best the data for optical conductivity and Raman 
scattering. The agreement with the data is quite good, 
not only in the positions of the maxima in optical con- 
ductivity and Raman response, but also in the shapes of 
(Ji(<jj) and R(uj). We argue that this good agreement is a 
strong argument in favor of a non-monotonic d x 2_ y 2 gap 
in electron-doped cuprates. 
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